Shining light through the Sun 
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It is shown that the Sun can become partially transparent to high energy photons in the presence 
of a pseudo-scalar. In particular, if the axion interpretation of the PVLAS result were true then 
up to 2% of GeV energy gamma rays might pass through the Sun, while an even stronger effect is 
expected for some axion parameters. We discuss the possibilities of observing this effect. Present 
data are limited to the observation of the solar occultation of 3C 279 by EGRET in 1991; 98% C.L. 
detection of a non-zero flux of gamma rays passing through the Sun is not yet conclusive. Since the 
same occultation happens every October, future experiments, e.g. GLAST, are expected to have 
better sensitivity. 
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Many contemporary models of particle physics predict 
the existence of light pseudoscalar bosons which are often 
refereed to as axion-like particles, or simply axions. 

The Lagrangian density of the photon-pseudoscalar 
system is given by 
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where F^ v is the electromagnetic stress tensor and F^ v — 
£fj,vp\F p \ its dual, a the pseudoscalar (axion) field, m the 
axion mass and M is its inverse coupling to the pho- 
ton field. The coupling of the photon and pseudoscalar 
fields in this way means that a photon has a finite prob- 
ability of mixing with its opposite polarisation and with 
the pseudoscalar in the presence of an external magnetic 
field [l|, l2| . One way of searching for such interactions is 
to pass photons through a magnetic field within which 
some conversion into pseudoscalars would be expected 
given a Lagrangian such as (JTJ) . The beam is then di- 
rected through an opaque medium, a wall for instance, 
and then passed though another magnetic field where any 
pseudoscalars which were created may convert back into 
photons. Experiments which aim to 'shine light through 
walls' in this way are more likely to succeed if one uses 
higher energy photons and stronger magnetic fields. In 
this letter we propose using the Sun as our magnetic field 
and as our wall, and gamma-ray sources from outside the 
solar system for our photons. This is feasible because the 
Sun does not appear to emit gamma-rays itself (at least 
in its quiet phase) 

Recently, the interest in axion-like particles has been 
re-ignited by the results of the PVLAS experiment which 
has detected a rotation in the plane of polarization of a 
laser passing through a magnetic field [4( ■ The level and 
nature of the effect is different from expected quantum 



electro-dynamics effects. It is claimed that the PVLAS 
result is compatible with the existence of such a pseu- 
doscalar with a mass of m ~ 1CV 3 eV and an inverse 
coupling M ~ 10 5 GeV. This is unexpected since the 
CAST |5| experiment has ruled out such a strongly cou- 
pled axion. Some authors have therefore suggested that 
there may be a composite pseudoscalar [6] with a form 
factor which suppresses its thermal production in the core 
of the sun while leaving its classical mixing properties 
unchanged. Since CAST looks for thermally produced 
axions, the compositeness is seen as a way of reconciling 
the results from CAST and PVLAS (see alsoQ). 

It is feasible that in a very few years, the axion inter- 
pretation of the PVLAS result may be tested using x-ray 
lasers and de-commissioned magnetic fields from particle 
accelerators @. Here we suggest that it might be tested 
on a similar time scale by searching for gamma rays from 
astrophysical objects when they are eclipsed by the Sun. 
(see also @, GjJ for related work) 

We are interested in spatial variations in the photon 
state vectors so we expand the photon field A(t, x) in 
components of fixed frequency A(x)e~ tut . If the mag- 
netic field changes on length scales much larger than 
the wavelength of the particles and the refractive in- 
dex | to — 1| <c 1 then one can expand Q the oper- 
ators in the equation of motion for the three species 
u) 2 + d 2 — > 2us(bj — id z ) and write down the Schrodinger 
equation 

id z -* = -{lo + M)^ ; V = (A x ,A y ,a) (2) 
where the matrix M. is given by 
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FIG. 1: Electron density and magnetic field near the surface 
of the Sun and into the convection zone. 
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FIG. 2: Axion parameter space. For the shaded region, GeV 
photon flux through the Sun detectable with GLAST at the so- 
lar occultations of 3C 279 is expected. The ellipse corresponds 
to PV 'LAS-motivated parameters Q/. The values of M above 
the thick line are disfavored by CAST f|/. 



the mixing parameters are the refractive index due to 
electrons in medium A p , the mass term for the pseu- 
doscalar A m and the off-diagonal which give the strength 
of the mixing Aj\/. The three parameters take values 
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where lu p — Ai:an e /m e is the plasma frequency, n e is 
the electron density, B is the magnetic field, m e is the 
electron mass, a is the fine-structure constant, u> is the 
photon's (axion's) energy. For constant magnetic field 
and electron density, the axion-photon conversion prob- 
ability is given by 
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where the oscillation length is 
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and z is the coordinate along the path. The resonant 
conversion takes place at w p ~ m, when the conversion 
probability is of order one provided the size of the res- 
onant region Az res is much larger than the oscillation 
length l osc . 

The solar model we choose is a combination of two 
models, one by Spruit [ll| for the convection zone (we go 
no deeper here) and one by Vernazza [12| for the chromo- 
sphere and photosphere. The electron density is shown 
in figure [TJ The magnetic fields inside the Sun are poorly 



known. We assume that they are in equipartition with 
the bulk motion of the plasma rather than the electron 
pressure, which gives us a clearly defined average value of 
the magnetic field at each radius within the Sun. How- 
ever, the magnetic field inside the Sun is turbulent, and 
it is necessary to include this variation in the magnetic 
field along the line of flight of the photon/pseudoscalar 
to obtain realistic mixing. We therefore assume random 
magnetic fields in the x and y direction generated by 
the Fourier transform of a power spectrum with random 
phases. For this power spectrum, we choose a flat spec- 
trum with a short distance cut-off at 100 km. The vari- 
ance of the magnetic field is normalised to be equal to 
the magnetic field at each radius. We see that the elec- 
tron density changes by 8 orders of magnitude in a layer 
of some 2000 km below the solar surface, thus providing 
conditions for the resonance for a wide region of the ax- 
ion parameter space. On the other hand, this region is 
close enough to the surface so one may expect that the 
probability for a photon, created in the resonant region, 
to escape from the Sun, is significant. 

For the energy range we consider, namely gamma rays 
with E > 100 MeV, the interaction between the photon 
and the medium is dominated by the pair production 
on the nuclear electric field well described by the Bcthc- 
Heitler cross section. 

To obtain indications of the regions in parameter space 
where the maximal mixing is allowed, we estimate the 
location and size of the resonant zone in the adiabatic 
approximation with the help of Eq. (j4]) and then im- 
pose two requirements: that the resonance takes place 
in the optically thin zone and that Az rcs > ^ osc . Quite 
straightforwardly, one obtains that the maximal mixing 
is allowed for relatively strongly coupled axions of various 
masses (see Fig. [2]). 

Inside the Sun, away from the resonance, the conver- 
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sion probability is very small. Despite this, the high den- 
sity and almost immediate absorption of any occasion- 
ally created photon may result in cumulative decrease 
of the axion flux. In the same adiabatic approximation, 
for the gamma-ray energies E < 10 GeV and for the 
paths not going very deep into the Sun (impact param- 
eters d > 0.7Rq), this cumulative effect does not re- 
duce the number of axions in the beam by more than 
10%. This means that for the parameters specified, a 
non-negligible part of the initial photon-axion mixture 
would travel through the optically thick regions in the 
form of axions and would convert back to photons at the 
solar surface. It is therefore possible that a gamma ray 
source would be able to shine through the Sun, although 
in practice only a fraction of photons would emerge. 

To be more precise and to go beyond the adiabatic 
approximation, we integrate numerically the equations 
of motion along a path of several thousand kilome- 
tres through the Sun. Instead of explicitly solving the 
Schrodinger equation, it is better to go to the interac- 
tion representation and separate out the mixing that we 
are interested in from the normal propagating oscillation. 
This can be done by defining a density matrix p = <3>*'I' 
with an evolution equation 0, [HI 



id z p = [M,p] - iVp 
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where here we have simultaneously introduced a damp- 
ing matrix T> designed to take into account interactions 
between the photon part of the wavefuction and the par- 
ticles inside the Sun. This damping matrix takes the 
form 
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where the inverse decay length T — n e a and we again 
use the Bethe-Heitler cross section for the a. We assume 
that the beam entering the Sun has a ratio of 2:1 pho- 
tomaxion for the following reason. Figure [3] is a Hillas 
plot of the magnetic field of some typical gamma-ray 
sources vs. their size (the transient nature of gamma ray 
bursts means that one could not predict when to point 
the detector at the Sun so we have not included them on 
this plot.) All of the sources are larger in size than the 
oscillation length at maximal mixing (l osc ~ M / B) cor- 
responding to their magnetic fields and the coupling of 
the PVLAS photon. A more detailed analysis would look 
at the coherence length of the magnetic field rather than 
the total extent of the object, but for the smaller objects 
with higher magnetic fields, one would expect these two 
length scales to be closer together. The right hand axis 
of figure [3] shows the energy necessary for maximum mix- 
ing to be achieved. The source we will be looking at in 
this paper is a quasi-stellar object, in other words an ac- 
tive galaxy, so we expect there to be strong mixing and 
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FIG. 3: Hillas plot of the magnetic field of objects vs. their 
size. The oscillation length at maximal mixing for the PVLAS 
coupling M — 10 J GeV is plotted as a function of magnetic 
field (diagonal line). The right hand axis shows the photon 
energy required for maximal mixing as a function of magnetic 
field assuming the PVLAS coupling and mass m = 10 -3 eV. 



many oscillation lengths in the source, and the high en- 
ergy photons created should therefore be equal mixtures 
of both photon polarisations and the axion. 

Figure[4]shows the actual mixing in action along a path 
through the Sun at a depth of 0.95 Rq. In figure [5] we 
plot the probability of a photon emerging as a function 
of the maximum depth along the path through the Sun 
for an axion mass of 10 -3 eV and different values of the 
coupling M. The results of the numerical solution are 
in rather good agreement with the analytical estimates 
presented above. In particular, the probability of gamma 
rays emerging from the Sun when a bright source passes 
behind it may be as high as 2% for GeV energy photons, 
if the interpretation of the PVLAS signal as an axion is 
valid, and even larger for the parameters shown in Fig. [2] 

Let us turn now to observational constraints on gamma 
rays passing through the Sun. The only source in the 
EGRET catalogue [l4j occulted by the Sun is the quasi- 
stellar object 3C 279. EGRET observed the solar occul- 
tation of 3C 279 in 1991, in the viewing period 11.0. The 
QSO was in moderate state and was firmly detected in 
gamma rays during that viewing period [14| . 

The source was screened by the Sun for 8 hours 34 
minutes on October 8, 1991 [15;]. The minimal impact 
parameter was d — 0.75Rq. 

We analysed the EGRET data following [16]. The 
apparent flux of photons with energy E >100 MeV 
from the location of 3C 279 during the occultation is 
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FIG. 4: The probability of a photon emerging when a GeV 
PVLAS axion enters the Sun with an impact parameter 0.95 
Rq. Upper curve: unrealistic case with photon cross sec- 
tion neglected. The two resonant regions can be clearly seen. 
Lower curve: the same but including photon interactions. 



a non-zero point-source flux from the location of 3C 279 
during the occultation cannot be confirmed or excluded 
at the present level of statistics. Since such a non-zero 
flux could be the signal of a new elementary particle, 
this calls for future astronomical observations with more 
precise instruments. In particular GLAST is predicted 
to have a sensitivity of 3 x 10 _8 cm~ 2 s _1 for the period 
of 8 hours [l3|. By looking at the same source when it 
is eclipsed by the sun, GLAST will therefore be able to 
rule out complete transparency to gamma rays. GLAST 
could also observe the partial transparency predicted in 
this paper, confirming the interpretation of the PVLAS 
data as a new particle. 
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FIG. 5: The probability of a GeV photon emerging for pseu- 
doscalar mass m = 10" 3 eV and different inverse coupling 
scales M. The dotted line bounds the maximal depths still 
transparent for GeV photons. 



obtained from the analysis of the rest of the same viewing 
period when the source is not behind the Sun which is 
(8.6 ± 0.5) • 10~ 7 cm" 2 s _1 . This latter result, obtained 
from our own analysis of the data, is in a very good agree- 
ment with the value quoted in the 3EG catalog [lj] for 
this period, (7.94 ± 0.75) • 10~ 7 cm~ 2 s -1 . We see that 
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